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Abstract. The harmonic sawtooth map w (x) of the unit interval onto itself is defined. It is 
shown that its fixed points {x : w (x) = x} are enumerated by the n-th derivatives of a Meijer- 
G function and Lerch transcendent, serving as exponential and ordinary generating functions 
respectively, and involving the golden ratio in their parameters. The appropriately scaled Mellin 
transform of w (x) is an analytic continuation of the Riemann zeta function f (s) valid V — Re(s) ^ 
N. The series expansion of the inverse scaling function associated to the Mellin transform of w (x) 
has coefficients enumerating the Large Schroder Numbers Sn, defined as the number of perfect 
matchings in a triangular grid of n squares and expressible as a hypergeometric function. A finite- 
sum approximation to ^ (s) denoted by ^u, (N; s) is examined and an associated function x s) is 
found which solves the reflection formula {N; 1 — s) = X ^) Cm s). The function x 
is singular at s = and the residue at this point changes sign from negative to positive between 
the values oi N = 176 and A'^ = 177. Some rather elegant graphs of the reflection functions 
X (A; s) are also provided. The Mellin and Laplace transforms of the individual component 
functions of the inflnite sums and their roots are compared. The Gauss map h {x) is recalled so 
that its fixed points and Mellin transform can be contrasted to those of w{x). The geometric 
counting function Nc^. (x) = ^^^+^ _ ij of the fractal string Cw associated to the lengths 
of the harmonic sawtooth map components {w„ {x)}'^_-^ happens to coincide with the counting 
function for the number of Pythagorean triangles of the form {(a, b,b + 1) : (6+1) ^ x}. The 
volume of the inner tubular neighborhood of the boundary of the map dCw with radius e is shown 

to have the particuarly simple closed-form Vc^ (e) = '^'^"^'^'^ 2^^)"^^'*^^ where v (e) = 2e 
Also, the Minkowski content of is shown to be Mc^ = 2 and the Minkowski dimension to be 
Dciu ~ 5 ^'^'^ thus not invertible. The geometric zeta function, which is the Mellin transform 
of the geometric counting function Nc^{x), is calculated and shown to have a rather unusual 
closed form involving a finite sum of Riemann zeta functions and binomioal coefficients. Some 
definitions from the theory of fractal strings and membranes are also recalled. 
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1. Unit Interval Mappings 
1.1. The n-th Harmonic Sawtooth Function Wn{x). 

1.1.1. Infinite Sum Decomposition. Let the n-th component ?i;„ (a;) G [0, l]Va; G [0, 1] of the harmonic 
sawtooth function w{x) 9, 2.3] be defined as 

(1) Wn{x) ^ n{xn + X - l)x {x,I^) 
where 

(2) x(^,/,f)= ^ ^^""^n 

{ otherwise 

is the characteristic function of the n-th harmonic interval ( |132|) . By setting n = \_x~^\ as in (|134p 
we get the unit interval mapping 

w{x) ^wi^-i^{x) 

= [x'^\ (x [x-ij + .T- l)x 

(3) =T,'^lWn{x) 

= J2n=l ^i^n + X-1)X (x, In) 

= [x'^\ {x [x-^\ + X- 1) 
= X [x^-^l + X \_x~^\ — \_x~^\ 

As can be seen in Figure [1] w{x) is discontinuous at a countably infinite set of points of Lebesgue 
measure zero 

H = {y: liui^^y- w{x) ^ lim^^2^+ w{x)} 
= {0,i:neZ} 

The left and right limits at the discontinuous points are 

lim^_^gH- w{x) = 1 
lim^_>gH+ w{x) = 



(4) 
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'I I — I — r^T — I — I — r — I — r 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



{w(x):x=0.1..1} 



Figure 1 . The Harmonic Sawtooth map 



1.1.2. The Fixed Points of w{x) as an Iterated Function System. The iterates of the map 
(6) [w{x),w{w{x)),w{w{w{x))),w{w{'w{w{x)))), . . .] = (x) , w"^ (x) , (x) , w'^ (x) , . . . 



have the form 



(7) 



w^{x) = a,. — (x — c)br 



where {a^, bj. E Wi : r E IN} is a pair of integer sequences and c S IR is some constant. The sequence 
of quotients ^ converges rather quickly to the fixed value 



(8) 



limr^oo = a; - c Vx e [0, 1] 
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The explicit equation(soinetinies called Schroder's equation) for the fixed points of is 

FixJI, = {x : Wn{x) = x} 

= {x : n{xn + x — l)x{x, I^) = x} 
= {x : n{xn + x ^ 1) — x} 

(9) ^hm.^o ^^j"^"' 

, -1 -1 -d) 6-1 



n-1 2-6 



71^ -\-n—l 

where G is the Meijer-G function and Fixt^(x) is the generation function 
Fix^(x) = Er=i FixJi^" 

(10) ~ Sn=l n^i+^-l 

= To ((5 - a/5)$ (a^, 1, 1 - 0) + (5 - V5)$ (x, 1, 0)) 

where $(z,a, t>) is the Lerch Transcendent (|159[) . and (p is the Golden Ratio, which is the ratio 
of two numbers having the property that the ratio of the sum to the larger equals the ratio of 
the larger to the smaller. [3 Ch.XX][351 1.7] [H p.50][lS] The number can be called the "most 
irrational number" because its continued fraction expansion, given by iterations of the Gauss map 
([75)) . converges more slowly than any other number. The constant cj) satisfies the simple identities 

(11) 



0-1 



(12) (t)"^ - <t) - 1 ^ 

An interesting fact is that the density of a motif in a certain noncommutative space described in 
[251 5.1] must necessarily be an element of the group K + 

1.2. Integrals Transforms of w(x). 

1.2.1. Dirichlet Polynomial Series and the Mellin Transform ofw{x). The Mellin transform of the 
harmonic saw map w{x), multiplied by 

(13) t{s) = s-l±l 

s — 1 

is an analytic continuation of the Riemann zeta function C(s)V — Re(s) ^ N. This form of the 
zeta function, denoted by Cwis), is the infinite sum of the Mellin transformations of the component 
functions. 

M[wn{x); X —^ s] = f"i Wn{x)x^~^dx 



/" Wnix^^)x~''~'^dx 



= /p n{xn + X — l)x{x, I^)x'' ^dx 
(14) = f" n{xn + X — l)x''~^dx 



1+1 

ein(^)(i-^)_„e° '°(?n^)_,e'""(l) 
n^n-\-s(n-\-l)^ — (n+l)^n 
n ~^ —n{n-\-l)~^ —sn~^ 
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There is a conjugate pair of inverse branches of t(s) found by solving 

T^\t) ={s:T{s)^t} 

(15) 

2 2 2 

where T^^{t) and TZ^{t) denote the positive and negative solutions respectively. The coefficients 
in the series expansions are integers enumerating the large Schroder numbers Sn which count the 
number of perfect matchings in a triangular grid of n squares, named after Ernst Schroder ( 1841- 
1902). [HI A006318][l p.340][Tn] 



(16) 



Sn — 2F1 



1 2 



12 



where pFq is a hypergeometric function. We have 



lim 



(17) 



(18) 



lim 



t-fO ■ 




n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 




1 


1 


2 


6 


22 


90 


394 


1806 


8558 


41586 


206098 


1037718 


5293446 


27297738 


142078746 



^ — n(n+l) ^ —s 



The residue at the singular point s = 1 of r(s) is 

Res(T(s)MK(x);x->s]) = Res (-1^^ 

= Res ;_"^''°+^"'° ) 



(19) 



= 2 f" Wn(x)Ax 

= 2 f " n{xn + X — l)da; 



The infinite sum of the Mellin transforms multiplied by t(s) analytically continues C(s)V(— Re(s)) ^ N 

C-(s) =C(s) 

= t{s)M[w{x);x s] 
= t(s) w{x)x'^^^dx 

= /o^ b^^J b^^J + a; - 1) x^-Mx 
= '^feTEr=iA^K(^);:^->s] 



(20) 



^ En=i /J^ n{xn + x- l)2;^-id:E 



_ Y^oo s±i {_!}^ 

— -^s-l \^ s(s+l) 

Eoo n(n+l)~^— n"'"~'^4-sn~ 
n— 1 s— 1 
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Unlike the Mellin transform of the Gauss map (|78p . which must be muhipHed by the factor s then 
subtracted from before it equals ({s), the "harmonic sawtooth continuation" Ctu(s) of the zeta 
function ({s) has the fortuitous property that it equals ({s) after multiplying M[w{x);x — s] by 
t(s). This property of w{x) allows us to put t(s) inside the sum to get an expression, denoted 
by Cu]{N; s), involving the difference of two Dirichlet polynomials, one of which is scaled by s. The 
substition oo — >■ A'^ is made in the infinite sum appearing the expression for Cto(s) to get 

(21) _ 1 f„ , /AT , tM-s 1 , „Y-^ „-s - 



~ {s~l){N+iy r(s) '> ^* 

with equality in the limit except at the negative integers 

(22) W->ooC«.(^;s) =C(s)V-s^N* 
The functions Cu, {N; s) have real roots at s = and s = — 1. That is 

(23) lim C^, {N; s) = hm Cw [N; s) = 

s— >0 s— > — 1 

The residue of {N; s) at s = 1 is given by 

Res(C^,(iV; s j) = Res (t{s) •£n=i M[wn{x); x ^ s]) 

S—1 S—1 \ / 

(24) = ELi Rff I T{s)M[wn{x);x s] J 



n 



2—m=l n'^+n 



N 
N+l 



Thus, as required 
(25) 



limAT^oo Res(C«, (N; s)) = limAr^oo E^Li WTr 



= 1 

The function t(s) has zeros at —1 and and a simple pole at s = 1 with residue 

(26) E'f(-(^)) =?if(^ST)=2 

The Mellin transform of r(s) has an interesting Laurent series, convergent on the unit disc, given 

by 

M[T{s)-S-^t] = T(s)s*-lds 

(27) =tn=rmn-2)e--^ y\t\<i 

sin(7rt) 
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The transformations M[wn{x);x — ^ s] have removable singularities at —1 and where the limits 
are given by 



lim^M[wn{x);x s] ^n^ln{^) +nln{^) 



^^^^ \imM[wn{xy,x ^ s] =l + nln(^) 

lims->o e-*=f["'"(^)'^^''l ={n + l)"7i-"e-i 

(29) lim,^_i e-A^[«'"(^);^^«l = n("'+")(n + l)(-"'-«)e" 

So, (IT4l) can be rewritten as 



n2ln(s±i) +nln(s±i) -n s = -1 
(30) M K(a;); a; ^ s] = <( 1 + nln{^) " s = 



otherwise 



Furthermore, we have the limits 



lim„^oo lim M [wn{x);x -> s] = lim„^oo lims^_i 

s— 1 



n — n 



(n+1) 



(31) =lim„^oon2ln(Ii±i)+nln(s±i)- 

_ 1 

2 

and 



lim„^oo lim M [wn {x);x s] = lim„^oo lim^^o - "''"^^^ 

(32) =lim„_ool + nln(^) 

= 

The sum of limits over n at s = — 1 is Euler's constant. [T5]|141 1.1] 

n— 1 n n 

= lims->iC(s) - ^ 

= 7 

0.577215664901533 . . . 

1.2.2. The Reflection Formula for C^, (iV; s). There is a reflection equation for the finite-sum ap- 
proximation Ctu(Af;s) which is similiar to the well-known formula (^{1 — s) = x{^)C{s) with 
X (s) = 2 (27r)"" cos {^) T (s). The solution to 

(34) C^iN;l-s)^xiN;s)Cn.{N;s) 

is given by the expression 



-n° + (Ti + l)°-'-n + n°-' 
(35) — ^ + + 



_ (s-1) X:"^! -rt° + (K-|-l)''~^w+K°~^-K°~^s 

which satisfies 

(36) x(^;i-s) = x(^;s)"' 
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The functions x indexed by N , have singularities at s = 0. Let 



= E^Li " (In (" + !)- In H) 

■L ( Y^-^ ln(n)n^ — ln(n+l)n'^ — ln(n) 

OV^} - 2^„=1 n(n+l) 

=iE,ti«(ln(" + lf -ln(nf) 

then the residue at the singular point s = is given by the expression 
(38) 

Res(x(iV; s)) = -Res(x(A^; s)'^) 



(37) 



l+7+*(n+2)-wTT+i:"=l 



a{N)-N 

N ln(n)„^-ln(,^ + l)„^-l„(„) " ('"(r( W + 1 ) ) - ^ E ^= i r. (in ( „ + 1 ) ^ - In ( „) ^ ) ) 
■>(" + !) (W-E"^i ti(ln(n + l)-ln(Ti)))(W + l) 



lim Res(x(A^;s)) = 1 

N—Kx> s=0 



(E"=i"(M"+l)-lnW))-iV 

which has the limit 
(39) 

We also have the residue of the reciprocal at s = 2 

2N 

(40) Res(Y(A^;s)"i) =7^7-77, ^ 

which vanishes as N tends to infinity 

(4f) lim Res(x(^;s)"i) = 

Af-foo s=2 

As can be seen in the figures below, the residue at s = changes sign from negative to positive 
between the values of iV = 176 and N = 177. 



^-2*(l,Ar+l)+2C(2) 
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Figure 2. 



|Res(x(iV;s)) : iV = 1...25o| 
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Figure 3. 



{Res{x{N;s))-^ -.N 

(. s=0 

For any positive integer N, we have the hmits 



= 1 



.250 



(42) 



The hne Re (s) — ^ has a constant modulus 



hnis^o X {N; s) 


= C!0 


hm.^o S^X (N; s) 


= oo 


lim^^i X(^;s) 


= 1 


hm^^i X {N; s) 


= 


\ims^2 X {N; s) 


= 




= 



(43) 



There is also the complex conjugate symmetry 

(44) xiN;x + iy) ^ x{N;x-iy) 

If s = n €E N* is a positive integer then x (^5 can be written as 
yfiV-n) - C..(jV;i^») 

Z^m = l Z^fc=l n \k-lj 

(45) ~ JV cos(^„)>I.(„-l.W-H) X 

^ ' (ii-i)(jv + i)" r(n) ^'^y'-J 

= N cos(,^,^)-I.(,^-l,W-^-l) ^, s 

(„-i)(jv + i)" r(n) -i-^yi-J 

The Bernoulli numbers pT] make an appearance since 

(46) xiN;2n)C^iN;2n) ^ B,^ (N + if ^^Il+R + , _ , 
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The denominator of x (-^5 ti^-s the limits 
(47) 



limAT^oo Cw {N; n) = ( [n) 
lim„^oo Cw (N; n) =1 



Another interesting formula gives the limit at s = 1 of the quotient of successive functions 



(48) 



lims=i 



X{N+I;s) ^ (jV+2)JV(JV+l-a(jV+l)) 

x(N;s) (N+lY{N-a{N)) 

_ (jV+2)Af(jV+l-E"=/»(ln("+l)~ln("))) 

~ {N+lf{N-Y."L n(ln(n+l)-ln(«))) 




Figure 4. {x (iV; s) : s = 1 . . . 2, = 1 . . . 25} 



Let 

(49) 



v{s) = x(oo;s) 



C(i-^) 

C(s) 



Then the residue at the even negative integers is 



(50) 



C(i-») 



n even 



Res(Ks))= 3lCWU^- 

nodd 



1.2.3. The Laplace Transforms L[wn{x)]x — >■ s\. The Laplace transform L[wn{x)]x — )■ s] and its 
roots are calculated to shed light on the behaviour roots of the Mellin transforms M[wn{x)\ x ^ s] 
but it is unclear whether this is accomplished. The Laplace transform (|135p of the n-th component 
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Figure 5. {x{N;s) : s ^ ^ . . .2,N ^ 1 . . . 100} 



function is given by 



(51) 



L[wn{x); X ^ s] = Wn{x)e ^*dx 

= n{xn + X — l)xix, In )G^^''dx 
= f \ n(xn + a; — X)e'^^Ax 



1+1 

n(n+l)e '^+1 — (n^+n+s)e 
71 



There is a removable singularity at s = which has the limit 
L[u;„(a;);a; -> 0] = lims^o i[w„(x); x - 



^ w(rt+l)e "+i -(rt +n+s)e 



^ ' = iRes (t(s) A/[w„(a;); a; ^ s]) 

Additionally, 



1 

1n{n+\) 



Eoo 1 
n=l 2n(n+l) 

(53) = lim„^oo \un^M[wn{x)]x ^ s] 



lim„^oo hnis^-i -- 
1 

2 
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The roots of L[wn{x); x — > s] are enumerated by 

(54) Pw^i'm) = {s : L['w„{x);x s] = 0} = -n{n + l)W{m,-e^''-) - n{n + 1) 

where W{m, x) is the Lambert W fmiction (11461) and rn G n e IN. It can be verified that 

where —e~^ is expressed as the continued fraction via its quotient sequence 

(56) - e-^ = [-1, 1, 1, 1, 2, 1, 1, 4, 1, 1, 6, 1, 1, 8, 1, 1, 10, 1, 1, 12, 1, 1, 14, 1, 1, 16, 1, 1, 18, 1, 1, 20, ...] 

The roots (m) satisfy a functional reflection equation with respect to m 



(57) PwA^) =Pw„i-m-'^) 



—n(n + l)W{m,—e ^) — n{n + 1) = —n(n + l)W{—m — 1, —e ^) — n(n + 1) 

where x — Re(a::) — Im(x) denotes complex conjugation. The quotients of the roots of consecutive 
transforms is 

/co^ _ n+l 

Thus 

(59) lim„^o. -lim™^ = l 

1.2.4. The Roots pt{ (m) of M[wnix); x s]. Define 

M^„(s, n) = T{s)M[wn{x); x ^ s] 

_ K(7t+l)~°-K^~°+STt~° 
~ S-1 

as in (1211) and its infintc number of inverse branches (which are currently lacking closed-form 
expression if such a thing is possible except when n = 1), where the branches are indexed by m 

M7j(z,n, m) = {s : T{s)M[wn{x); x s] ^ z} 

(61) ^ f „ . n(ri+l)-°-n^-°+sn-° ^ ^\ 

then we see that the first function where n = 1, M^J(z, l,m), has the closed-form 
where is the Lambert W function (|146l) . It is verified that 



1-2 i"<2) — ^ ' :„,o>^ 



+ln(2) 



(63) M_(M-J(^,l,m),l,m) ^ ^y^,^ , ^ 

hT(2) 

Furthermore, let (m) denote the m-th root of M[wn{x); x ^ s] 

, . ptiim) = {s: M[wnix);x ^ s]=0} 

^ ' = Af-HO, n,m) 

which satisfies 

(65) Im {ptl (m)) > Im (^pt{_^ (m)) 
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Figure 6. The Roots /0^„(to) of {L[w„(x); x ^> s] : n £ 1 .. .9} 

(66) lini„^±oo Re(p^^„ (m)) =0 
Thus, 

(67) hm„,^±oo arg (p*^^ (m)) = f 

1.2.5. Quotients and Differences of (m). Let 

(68) Ap^Jm) ^P^Sm + l)-p^JJm) 

be the forward difference of consecutive roots of M[wn{x)] a: — >■ s]. The hmiting difference between 
consecutive roots is the countably infinite set of solutions to the equation rtt + (n + l)i = given 
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by 



(69) 



Ap^^ (±00) = lim^^ioo Ap^[ (m) 

= {s : n5 + (n + 1)5 = o} 
= lim„^±oo {p^Jm + 1) - p^Jm)) 

= limm^-ibcx) j^ ^""^" ^"* 

_ 27ri 

~ M[x(a;,/^);a;->0] 



Let denote the limit 



ln(n+l)-lii(n) 



0°° -lim ^ P-nC'") 

(±00) 









)-In(,i-l) 


l„(n + l)-ln(,.) 


h,(„,- 


l)-lu(n-2) 


, ln(,i 


)-ln(„-2) 


In(ri)- 


-ln(ri-l) ^ 



(70) ApM_J±oo) 

^ ' _ M[x(.,(^.^));.^0] 

M[x(a:,(i,^));a;^0] 
_ ln(Tt)-ln(Tt-l) 
~ ln(n+l)-ln(n) 

then we also have the limit of the limits Q^m as n — > 00 given by 

lim„^.±oo QpM = hm„_^±c 

(71) = lim„_>±c 

= lim„^.±c 
= 1 

— p^ (m+1) M / N_ M / 

The limiting quotients ^ = e^^n'-^'' ^^^n*-^"^ ^ as m ^ oo are given by 

(72) 

lim„^^ ef"n(m)-P^n(m+i) = 1 - 2sin - 2icos { ,^^^^,]_,^^^^ ) sin 

27ri 

= e ln(n+l)-ln(») 

where we have 

I lim e'"»n(™)-''»n(™+i)| -lim |e''"n('")-''»n('"+i) I 

I i'i"m— >CX) en- n- I — 111I1to_).oo |e n- n- | 

(73) _ -\/^-iiii3f5i(i7 

= 1 

and 

(74) lim„^^ lim„^^ eP^„(™)-C(™+i) = -1 



INTEGRAL TRANSFORMS OF THE HARMONIC SAWTOOTH MAP, THE RIEMANN ZETA FUNCTION, FRACTAL STRINGS, AND A FINITE 




Figure 7. {pi^^(m) : m = 1 . . . 5} 



1.2.6. The Laplace Transform L [{s — l)M[w„(x); x -> s] ; s — > t\. The Laplace transform of (s — 
l)M[wn{x); X — > s] defined by 

L[{s — l)M[wn{x);x — > s] ; s t] = L[n{n + 1)^^ + sn"" — n^""; s — > 

t+ln( )t+ln(»+l)+7i ln(ra)^-n ln(n) ln(n+l) 

^ (ln(n)+t)2(ln(n+l)+t) 

has poles at — ln(n) and — ln(n + 1) with residues 

Res ( LUs — l)M[wn(x); X ^ s]; s ^ t]] = —n 

t=- ln(n) V / 

(76) / X 

Res [L\{s - l)M\wn(x);x ^ s\]s ^ t]] 

t=-ln(n+l) V V 



1.3. The Gauss Map h{x). 

1.3.1. Continued Fractions. The Gauss map h{x), also known as the Gauss function or Gauss trans- 
formation, maps unit intervals onto unit intervals and by iteration gives the continued fraction 
expansion of a real number [38l A.1.7][35l I.l][l21 X] The rt-th component function hn{x) of the 
map h{x) is given by 

(77) Kix) ^l^xi^J^) 
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The infinite sum of tfie component fmictions is the Gauss map 



(78) 




{h(x):x=0.1..1} 



Figure 8. The Gauss Map 



The fixed points of h{x) are the (positive) solutions to the equation hn{x) — x given by 

Fix^ = {x : K{x) = x} 

(^9) ^{x:l^ = x] 

■\/k^+4 n 

~ 2 1 
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1.3.2. The Mellin Transform of h{x). The Mellin transform (|138p of the Gauss map h{x) over the 
unit interval, scaled by s then subtracted from is an analytic continuation of ({s), denoted by 
Ch{s), valid for all— Re(s)) ^ N. The transfer operator and thermodynamic aspects of the Gauss 
map are discussed in [42] [41] [40] [39] [36] . The Mellin transform of the n-th component function 
Wn (x) is given by 

M[/i„(x); x — > s] = j}^ hn{x)x'^~^dx 

^loJ^X {x,I,'^)x^-'dx 
(80) =/i_ [x-ij)a:^-idx 

= f\ l^x'-^dx 

_ Tt(n+1) "+s{n+l) °-n-^~° 



which provides an analytic continuation Ch{s) — C('5)V(— Re(s)) ^ N 

Chis) = - sM[h{x);x s] 
= 7^-sfoh{x)x'''^dx 
(81) =i3T-s/o (2;"'- b"'J)^'"'dx 

= j^-s X:r=i M[K (x) ; X s] 



1.4. The Harmonic Sawtooth Map w(x) as an Ordinary Fractal String. 



1.4.1. Definition and Length. Let 

(82) X,f -{^^^^) 

be the n-th harmonic interval, then {w{x) G C^^ : x G fl} is the piecewise monotone mapping of 
the unit interval onto itself. The fractal string £^ associated with w{x) is the set of connected 
component functions Wn{x) C wix) where each w„(x) maps/^ onto (0,1) and vanishes when 
X ^ . Thus, the disjoint union of the connected components of Cw is the infinite sum w(x) = 
^^j^ Wn{x) where only 1 of the Wn{x) is nonzero for each x, thus w{x) maps entire unit interval onto 
itself uniqely except except for the points of discontinuity on the boundary dCw = {0, — : n £ N*} 
where a choice is to be made between and 1 depending on the direction in which the limit is 
approached. Let 

(83) ujn{x) — n{xn + X — l)x{x, I^) 
where xi^jln) is the n-th harmonic interval indicator (jl27p 

(84) x(^,/,f) = (^1^) - ^ (^) 
The substitution n — >■ [ij can be made in (|132p where it is seen that 
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and so making the same substitution in (155)) gives 



(86) 



= Er=l ^i^n + ^)X{X, In) 

= [x-^\ {x [x-^\ +X-1) 




Figure 9. The Harmonic Sawtooth Map 



The intervals 1'^ will be defined such that £wn = |/™| = |t«ri(x)|. Let 
(87) 



f]n = JjH x{n + l)nAx 

_ (^T±i) 



2 

2n+l 
2n(n+l) 
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be the midpoint of then 



_ ( 4n+l 4n+3 
4n(n+l) ' 4n(n+l) 



SO that 



= |w„(a;)| 

= /n^?T.(a;n + x- l)x{x,I^)Ax 
= f "1 w(a;)dx 
(89) = |/H 



/o'x(a:,/;r)dx 

4w+3 _ in+1 
An{n+l) in(n+l) 

2n(n+l) 




The total length of is 



(90) 



Eoo 1 
r. 

_ 1 
~ 2 



= 1 2n(n+l) 
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1.4.2. Geometry and Volume of the Inner Tubular Neighborhood. The geometric counting function 
(IT091) of is 



(91) 



Nc^ix) = #{n ;? 1 : £w-^ ^ x} 

= 4{n ^ 1 : 2(n + l)n ^ a;} 



V2a: + 1 1 
2 2 



which is used to calculate the limiting constant (lllip appearing in the equation for the 
Minkowski content 



(92) 



= Um^;- 

2 



2 2. 



The function Nc„{x) happens to coincide with [T71 A095861], which is the number of primitive 
Pythagorean triangles of the form {(a, 6, & + 1) : (6 + 1) n}. ^ 171-176] [371 10.1] [H 11.2-11.5] 
Let 



(93) 



v{e) — min(j : Iwj < 2e) 



2e 
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which is the floor of the solution to the inverse length equation 

= {n : tw,,-i = 2e} = [n : = 2e] 



(94) 



tw 



= £ 



Then the volume of the inner tubular neighborhood of dCw with radius e pOSp is 






(95) 



1 

^ v{e) 2(n+l)n 



2e 




1 


2 


2 


4£V 


'ef-iev 


e) + l 



1 

2^ 



since 
(96) 

and by defintion we have 
(97) 



Eoo 1 
n=m 2n(n+l) 



1 

2m 



\iTn,^^^Vc^{e) =0 
linie^oo Vc„(£) = |£^| = i 

Thus, using ^ and ([25]), the Minkowski content piU)) of is 



: lim„ 



(98) 



_ ^^2^-^ 





1 








2 2 




2e 





= 2 



1.4.3. The Geometric Zeta Function C,c^{s). The geometric zeta function (jll2l) of w{x) is the 
Dirichlet series of the lengths ((M)) and also an integral over the geometric length counting 
function ^ Nc^{x) 



(99) 



En=l (2n(n+l)) 



V2^+T 
2 



idx 

.T^^'^-^da; 
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Figure 12. Geometric Counting Function Nc,^{x) of w{x) 



The residue pl5)) of (s) at D^^ is 

Res(C£(s)) -lim^^^+ (s^DcXds) 



(100) 



= lim 

s 

— lim 

s 

— lim 

s 

= 



.1+ (■'^-5)C£W 

(.s-^)Er=l2-^(r. + l 

4+ («- 5)^/0 



V2a;+1 
2 
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The values of (c^ {n) at positive integer values n G N* are given explicitly by a rather unwieldy sum 
of binomial coefficients and the Riemann zeta function C{n) at even integer values. First, define 



(101) 



bn 
dn 



^ («-l)(l-(-l)"+i) 

= (-l)"(n-l) 

_ (-1)" 
2 



(-1) 
4 



then 
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Figure 13. {v{e) = mm{n : £wn < 2e) : e = ■ ■ ■ \} 



The terms of (n) from n 
( +\ 



35 

M 

'Mi 
1^ 



(^r(^ b„ 2(-ir( 



27n + c„-d„ + i 



1 to 10 are shown below in Table [TJ 



)C(<i„+2n-|-2m-c„) 



(^35 

4§j§9 
512 



+ U(2) 






-|C(2) 






+|C(2) 
-|C(2) 


+ K(4) 




-|C(4) 
+iC(4) 
-§C(4) 
+ %C(4) 
-#C(4) 
+ #C(4) 




+ iC(2) 


+ ^C(6) 
-SC(6) 


+ f C(2) 


-lie (6) 
+ SC(6) 


-^C(2) 
+ #C(2) 



-t|8C(8) 
-4C(8) 

Table 1. {US^) = S row„:n=1..10} 
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Figure 14. Volume of the inner tubular neighborhood of dCw with radius e {Vc^ (e) : e = . . . |} 

2. Fractal Strings and Dynamical Zeta Functions 

2.1. Fractal Strings. A a fractal string C is defined as a nonempty bounded open subset of the 
real line £ C R consisting of a countable disjoint union of open intervals Ij 

oo 

(103) = U 

The length of the j-th interval Ij is denoted by 

(104) £,^\I,\ 

where | • | is the 1-dimensional Lebesgue measure. The lengths must form a nonnegative monot- 
ically decreasing sequence and the total length must be finite, that is 



(105) ^J=i ^ 



l^|i=E,°:i^.<oo 
^1 ^ 4 ^ . . . ^ ^ ij+i ^ • • • ^ 



The case when = for any j will be excluded here since ij is a finite sequence. The fractal string 
is defined completely by its sequence of lengths so it can be denoted 

(106) C 
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0.02 0.04 0.06 0.08 0.10 0.12 



Figure 15, {li^ : e = O. , .i} and ^^^^ = \/2 

The boundary of L in R, denoted by dL C fi, is a totally disconnected bounded perfect subset 
which can be represented as a string of finite length, and generally any compact subset of R also 
has this property. The boundary dZ is said to be perfect since it is closed and each of its points 
is a limit point. Since the Cantor-Bendixon lemma states that there exists a perfect set P C dH, 
such that dL — P is a most countable, we can define £ as the complenent of dC in its closed convex 
hull. The connected components of the bounded open set are the intervals Ij. [551 1-2] [311 

2.2 Exi7]^, 3.i][5S][ii][is][in][ii][n][i3[2n] 



2.1.1. The Minkowski Dimension Dc and Content Mc- The Minkowski dimension Dc G [0,1], 
also known as the box dimension, is maximum value of V{e) 

(107) Dc -inf{a^0:V^(e) = O(ei-")ase^0+} = C£(l) = Er=i^j 
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where V{e) is the volume of the inner tubular neighborhoods of dC with radius e 

Vc{e) ^\xeC:d{x,dC)<e\ 

ij^2e ij<2e 

= y 2e+y 

(108) , , 

j 

and A'£ (x) is the geometric counting function which is the number of components with their recip- 
rocal length being less than or equal to x. 

Ncix) =#{j^l:£-' ^x} 

(109) _^^J'^- / 



The Minkowski content of C is then defined as 

Mc =lim,^o+^^ 



(110) 



where Cc is the constant 



1-D£ 

Dc{l-Dc) 



(111) Cc= lim.^oo^ 

If A4c G (0, oo) exists then C is said to be Minkowski measurable which necessarily means that the 
geometry of £ does not oscillate and vice versa. [211 1] [3][2l][26l 6.2] 

2.1.2. The Geometric Zeta Function Cc{s). The geometric Zeta function Cc{s) of C is the Dirichlet 
series 

^ ' ^sJ^NcJx)x-^-'dx 

which is holomorphic for Re(s) > Dc- If C is Minkowski measurable then < Dc < 1 is the simple 
unique pole of Cc{s) on the vertical line Re(s) = Dc- Assuming Cc{s) has a meromorphic extension 
to a neighboorhood of Dc then Cc{s) has a simple pole at Cc{Dc) if 

(113) Ncis)^0{s^^)ass'^oo 
or if the volume of the tubular neighborhoods satisfies 

(114) V^£(£) = 0(ei--°^)ase^O+ 

It can be possible that the residue of Cc{s) at s = Dc is positive and finite 

(115) 0< lim {s - DcKcis) < oo 

even if Nc{s) is not of order s^'^ as s — >■ cx) and Vc{s) is not of order £^~^^, however this does not 
contradict the Minkowski measurability of C. 
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2.1.3. Complex Dimensions, Screens and Windows. The set of visible complex dimensions of £, 
denoted by VciW), is a discrete subset of (D consisting of the poles of {(^(s) ■ s e W}. 

(116) Vc{W) ^{weW ■.Cc{w)\saY>o\e} 

When W is the entire complex plane then the set I?£((D) = Vc is simply called the set of complex 
dimensions of C. The presence of oscillations in V{e) implies the presence of imaginary complex 
dimensions with Re(-) = Dc and vice versa. More generally, the complex dimensions of a fractal 
string C describe its geometric and spectral oscillations. 

2.1 A. Frequencies of Fractal Strings and Spectral Zeta Functions. The eigenvalues An of the Dirich- 
let Laplacian Am(x) — —-^u{x) on a bounded open set £7 C IR correspond to the normalized fre- 
quencies /„ — of a fractal string. The frequencies of the unit interval are the natural numbers 
n S N* and the frequencies of an interval of length i are ni~^ . The frequencies of Z are the numbers 

(117) A,, = fc^7Vfc,j€N* 

The spectral counting function N^n (x) counts the frequencies of £ with multiplicity 

(118) =i:r^i^-(f) 

The spectral zeta function Ct'£(s) of C is connected to the Riemann zeta function (??) by 

(119) =c(s)Er=i^j 

= C(5)C£(5) 

[m i.i][isi 1.2.1] 

2.1.5. Generalized Fractal Strings and Dirichlet Integrals. A generalized fractal string is a positive 
or complex or local measure r]{x) on (0, oo) such that 

(120) / ■q{x)dx = 



for some 2:0 > 0. A local positive measure is a standard positive Borel measure 77(7) on (0, cxd) 
where J is the set of all bounded subintervals of (0, cxd) in which case r]{x) = |??(a;)|. More generally, 
a meausre r]{x) is a local complex measure if ri{A) is well-defined for any subset A C [a, b] where 
[a, b] C [0, 00] is a bounded subset of the positive half-line (0, 00) and the restriction of ry to the 
Borel subsets of [a, b] is a complex measure on [a, b] in the traditional sense. The geometric counting 
function of r]{x) is defined as 

(121) N^ix) = J^vix)dx 

The dimension is the abscissa of conergence of the Dirichlet integral 



(122) C|,|(^) = / ^'''\vix)\dx 

Jo 

In other terms, it is the smallest real positive a such that the improper Riemann-Lebesgue converges 
to a finite value. The geometric zeta function is defined as the Mellin transform 



(123) Cnis) = / x-''r]{x)dx 

where Re(s) > Djj. 
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2.2. Fractal Membranes and Spectral Partitions. 

2.2.1. Complex Dimensions of Dynamical Zeta Functions. The fractal membrane Ti associated 
with C is the adehc product 

oo 

(124) n^JlT, 

where each Tj is an interval Ij of length \og{i~^)~^ . To each Tj is associated a Hilbert space 
Hj = L^{Ij) of square integrable functions on Ij. The spectral partition function Zc{s) of £ is a 
Euler product expansion which has no zeros or poles in Re(s) > Dm{^)- 

where Dm{C) is the Minkowski dimension of C and Zcj{s) = is the j-th Euler factor, the 

j 

partition function of the j-th component of the fractal membrane. 23, 3.2.2] 

2.2.2. Dynamical Zeta Functions of Fractal Membranes. The dynamical zeta function of a fractal 
membrane C is the negative of the logarithmic derivative of the Zeta function associated with C 

Zeis) ^-4,HCc{s)) 

(126) _ ^U{s) 

3. Special Functions, Definitions, and Conventions 
3.1. Special Functions. 

3.1.1. The Interval Indicator (Characteristic) Function x{x, I) . The (left-open, right-closed) inter- 
val indicator function is xi^i ^) where / = (a, 6] 

1 X e I 
X ^ I 

(127) _ / 1 a<xs^b 

otherwise 

^ 0ix - a) - e{x - a)0{x - b) 

and 9 is the Heaviside unit step function, the derivative of which is the Dirac delta function S 

jS{x)dx ^0{x) 

(128) _ r x<0 

" \ 1 x^O 

The discontinous point of 0{x) has the limiting values 

^29) limj^^o- ^ 

lim2,^o+ d{x) 1 

thus the values of x(a;, (a, b)) on the boundary can be chosen according to which side the limit is 
regarded as being approached from. 

lim^^a- X{x,{a,b]) =0 

^3QN lima;^a+ Xix, (a, b]) = 1 - e{a - b) 

lim^^b- x{x, (a, b]) = 6{b - a) 

lim^^fc+ X(a:, (a,6]) =0 
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3.1.2. "Harmonic" Intervals. Let the n-th harmonic (left-open, right-closed) interval be defined as 



(131) 

then its characteristic function is 
(132) 



H 



1 1_ 

n+1 ' n 



-1) - 0(^-1) 
n+1 / V n/ 

9 {x — n) 



n+1 ^ 

1 -^<Xi^ 
n+1 ^ 

otherwise 



As can be seen 
(133) 



I \°° I" 

Un=l ^n 
X/n=l X(^5 ^n ) 



UZi{^i^^] =(o,i] 



The substitution n — > [ij can be made in (|132p where it is seen that 



(134) 



a; [x """J+x— 1 



1 Va;e[-l,-fl] 



3.1.3. T/ie Laplace Transform L'^[f{x);x — > s]. The Laplace transform 1.5] is defined as 
(135) Li[f{x);x^s] ^ Jlf(x)e-^dx 

where the unilateral Laplace transform is over the interval (a, b) — (0, oo) and the bilateral transform 
is over (a, &) = (—00,00). When (a, 6) is not specified, it is assumed to range over the support of 
f{x) if the support is an interval. If the support of f{x) is not an interval then (a, b) must be 
specified. Applying L to the interval indicator function (jl27p gives 



I-aixix, {a,b));x 



(136) 



= j\0{x - a) - e{x - b)0{x ~ a))e-^Ma; 



-s(!, + a) 



The limit at the singular point s = is 



(137) 



lim^^o Ll[x{x, (a, 6)); x ^ s] = limg^o 



(e'"-e''°)e 



3.1.4. The Mellin Transform M^[f{x);x-+ s] 
defined as 

(138) 



The Mellin transform [3U 3.2] [33 n.l0.8][5l 3.6] is 



M'Mixy,x^s] ^ jl;f{x)x^-^dx 

where the standard Mellin transform is over the interval (a, b) = (0, 00). Again, as with the notation 
for Laplace transform, the integral is over the support of f{x) if the support is an interval and (a, b) 
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is not specified, otherwise (a, b) must be specified. Applying M to the interval indicator function 
([T27| gives 



(139) 



M^,[x{x, (a, 6)); x ^ s] = J° x{x, {a, b))x'^-^dx 



J {e{x - a) - 9{x - b)9{x - a))x^-idx 



The limit at the singular point s = is 

[Xix, (a, b));x^O] = M [x{x, (a, 6)); a; ^ 0] 

= lim.^o M^[x{x, (a, 6)); x ^ s] 
= lim.^o ^ 



(140) 



= ln(6) - ln(a) 

The Mellin transform has several identities [31] 3.1.2], including but not limited to 

M[f{ax);x s] = a-^M[f{xy, x ^ s] 



(141) 



M[x°'f{xy,x ^ s\ 
M[f{x°');x~^ s] 
M[f {x-°');x ^ s] 

M[a:"/(a;^);x ^ s] 

M[a;"/(a;-^);x ^ s] 
A/[ln(x)"/(a;);a; ^ .s] 



M[f{x);x -> s + a] 
lM[fix);x^i] 
iM[/(x);x. 



/(.T);.T^^ 



-g;Af[/(x);x^s] 



where a > 0, /i > 0, and n € N. The Mellin transform of the harmonic interval indicator function 

una) is 



M 



(142) 



[X (x, ) ; X ^ 0] = X (x-, (;^, i) ) x^-^dx 



which has the limit 



(143) 



s 

M [x {x, I^) ; X ^ 0] = lini.^o M [x {x, I^)]x^ s] 
= ln(n + 1) — ln(n) 



The Mellin and bilaterial Laplace transforms are related by the change of variables x 
resulting in the identity [ST, 3.1.1] 



My) 



(144) 



Mo-[/(-ln(x));x ^ s] - LtZ[f {v): V ^ s] 
/(- ln(x))x-Mx = f{y)e'y^dy 



3.1.5. The Lambert W Function W{k,x). The Lambert W function [S][33] is the inverse of xe^ 
given by 

W{z) = {x : xe^ ~ z} 
= W{0,z) 



(145) 



l + (ln(z)--l)exp(^/„°°^ln 

EC30 
fe=l 



dx 



k\ 
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where W{a, z)\/a ^TL,z ^ {0, -e ^} is 

(146) W{a,z) = 1 + (2z.a + ln(.) - 1) exp ^ In [ Z% 
A generaliztion of ()145p is solved by 

(147) {x:xh- = z] 

The W function satisifes several identities 

M^(z)e^(^) = z 

w\zh\{z)) =ln(z) Vz < 1 

\W{z)\ =W{\z\) 

(^^^^ \n{W{n,z)) =\n{z)--Win,z) + 2iTTn 

_ \n{z) 



Wl^^'-^j =-ln(z) Vze[0,e] 

W{-ln{z)) _ z'^' 

-\n{z) 2 

where n G Z. Some special values are 

VK(-l,-e-i) =-1 

W{-e-^) = -1 

Wle) = 1 

T4^(0) = 

(149) W(oo) = oo 

W^( — OO) = oo + ZTT 

W{-^) =f 
W{-\n{V2)) =-ln(2) 
M^(-l,-ln(V2)) =-21n(2) 

We also have the hmit 

(150) lini,^±^i5^ 
and differential 

^ ' dz^^ y^-^Jy^)) - f(z)(l+W(a,fiz))) 

as well as the obvious integral 

(152) /^it^(_M£))dx ^/„^-ln(x)d2; = l 

Let us define, for the sake of brevity, the function 

W^i„(^) =1^(-1,-^) 

^ 1 + (>" (-^) - ' - ^"-) - r A ... (SSW)''^; 

Then we have the limits 

/-,,.-, lim:r^-oo Wln(x) =0 

lim^r^+oo Win(a;) = -oo 
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X = - e ..2 e 



■g{(-ln(x)) ^[-Hx)) |-ln(x)| S«(^inW) 



Figure 16. W 



(-^) = - and W^i„(a;) = W/ (-1, -IhM^ 



and 



(155) 



Im(W^i„(a:)) 



— TT — cx) < a; < 
... < a; ^ 1 
1 < a; < oo 



(156) Win{x) =-\n{x) yx^[0,e 

The root of Re {Win{x)) is given by 

{a::Re(m„(a:))=0} 



(157) 



{-x: {x^)~ =e3-} 
^ttW (Itt) 
0.27441063190284810044 . 
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where the imaginary part of the value at the root of the real part of IVin(z) is 



(158) =w(-l.Sf) 

9i il.712388980384689857. 



3.1.6. The Lerch Transcendent $(z,a, w). The Lerch Transcendent '14', 1.11] is defined by 



(159) Hz,a,v) V{|z|<l}or{|z| = landRe(a)>l} 



The Riemann zeta function is the special case 



(160) C{s) =ms,i)^En=ori' 



3.2. Applications of w(x). 



3.2.1. Expansion ofj. Consider Euler's constant 7 = 0.577215664901533 . . . ([33 
(161) u;"(7) =an-fo„7 



whereupon iteration we see that 
(162) 









1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


i 




• 


1 


48 


290 


581 


1163 


2327 


13964 


7492468716 


14984937433 


1078915495184 








2 


84 


504 


1008 


2016 


4032 


24192 


12980362752 


25960725504 


1869172236288 
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3.3. Conventions and Symbols. Many of these symbols are from [331 p491]. 



(163) 
i 

M {x : — oo < a; < oo} 

M {x : — oo ^ a; ^ oo} 

]R+ {a; : ^ a; < oo} 

M'^ {xi . . .Xd '■ — oo < Xi < oo} 

C {x + iy : x,y eM} 

% {...,-2,-1,0,1,2,...} 

N {0,1,2,3,....} 

W {1,2,3,....} 

H {o,i:neS} 

!{x) = 0{9{x)) fg<(^ 

/(a;) = o{g{x)) lim.^oo fg} = 

/(a;)xg(x) {a «C : {a, 6} > o} 

#A numbers of elements in the finite set A 

\A\d (i-dimensionalLebesguemeasure(volume) of A C R'^ 

d[x,A) {min(|a; — y|) : y G A} Euclidean distance between a; and the nearest point of A 

exp(a;) exponential e"" = I]^o C" 

Res(/(a;)) complex residue of /(a;) at a; = y 

\x\ floor, the greatest integer ^ x 

{a;} a; — [xj , the fractional part of x 

X complex conjugate, Re(x) — Im(a;) 

Fix J n-th fixed point of the map / (a;) , n-th solution to f{x) — x 

Pk fc-th prime number 
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